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ABSTRACT 

High-strain piezoelectric materials are often ceramics with a complicated constitution. In particular, PZT is 
used with compositions near to a so-called morphotropic phase boundary, where not only different variants of 
the same phase (domains), but different phases may coexist. Micro- mechanical models for ferroelectric ceramics 
would be much more realistic, if these effects could be incorporated. 

In this paper, we consider the conditions of mechanical and electrical compatibility of ferroelectric domain 
structures. We are able to address the question of coexistence of different crystallographic phases within the 
very same crystallite. In general, the spontaneous strain and spontaneous polarization of different phases are not 
compatible. The numerical analysis of the derived relationships are susceptible to the crystallographic description 
of the phases in question. In this presentation, a simple analysis and analytical, composition dependent fit of 
strain and polarization of PZT at room temperature for available data are used. 

The outlined approach can be used to model the overall behavior of multi- variant and multi-phase crystallites 
with certain, simplified geometrical arrangements of the constituents. 
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1. INTRODUCTION 

High-strain piezoelectric materials are often materials with a complicated constitution. The most interesting 
materials, PZT (PbZr x Tii_ x 3 ), (PbMg 1/3 Nb 2/3 3 )i_ x - (PbTi0 3 ) x (PMN-PT) and PbZn 1/3 Nb 2/3 3 )i_ x - 
(PbTi0 3 ) x (PZN-PT) have in common, that they are solid solutions which exhibit a so called morphotropic 
phase boundary. These materials show extreme material properties, namely high piezoelectric coefficients and 
are therefore often used at compositions close to the morphotropic phase boundary (MPB). 

The MPB in the phase diagram (see fig. 1) separates phase regions of two different ferroelectric phases, in 
the case of PZT the tetragonal and the rhombohedral phase. The MPB was also considered to be a small region 
of coexistence of different phases. In regard of the exceptional properties of these materials, it is often discussed, 
that a phase transition of one phase to the other would be possible under applied electric fields, for instance. This 
would give the material a larger degree of freedom to respond to applied loadings and would explain exceptional 
electromechanical coupling phenomena. 

Less than ten years ago, a closer investigations of PZT at compositions at the MPB 3-5 have revealed, that 
there exists a small region, where a monoclinic phase can be observed. The existence of the monoclinic phase 
is more pronounced at low temperatures, but it is also observable at room temperature, mostly in coexistence 
with the tetragonal phase. For a review of the crystal structural studies see Ref. 1. It can also be deduced, that 
the monoclinic phase a intermediate phase between the tetragonal and the rhombohedral phase and that both, 
tetragonal space group P4mm and rhombohedral space group R3m are subgroups of the monoclinic space group 
Cm (see fig. 1). This opens the possibility that a tetragonal to rhombohedral phase transition could take place 
via a monoclinic intermediate state. 

More recently, new experiments brought new arguments into the discussion about the nature of the mor- 
photropic phase boundary. Glazer et al. 6 have concluded from local electron diffractions, that the local structure 
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Figure 1. left: Schematic phase diagram of PZT, after Ref.. 1 Besides the traditional cubic, tetragonal and rhombohedral 
phases, the phase diagram shows schematically also the recently discovered intermediate phases Cm, Pm and Pc. An 
additional low temperature phase (Cc as an intermediate phase between R3c and Cm) at the morphotropic phase boundary 
is not shown, right: Space group relationships among tilt-free space groups of ABO3 perovskites. Lines joining two space 
groups indicate a group-subgroup relationship, where solid lines stand for possible 2nd. order (continuous) transformations 
and dashed lines indicate first order transformations. After Howard and Stokes, Ref. 2 



of PZT is always monoclinic, and the average tetragonal and rhombohedral structures are a result of short-range 
to long-range ordered states. These new viewpoints may stimulate new constitutive models of ferroelectrics with 
compositions close to a MPB, but so far, a constitutional model, which deals with changes of the structural 
order is not available. Therefore, micromechanical models will still be a starting point to discuss the behavior 
of the material. This is also justified by the clearly visible ferroelectric domains which represent ordered states 
at a mesoscopic level. For that reasons, the purpose of this paper is to examine the compatibility of tetragonal, 
rhombohedral and monoclinic domains, depending on the chemical composition of PZT across the MPB. The 
starting point of our study is a homogenous crystal which is in a virtual cubic state (state of reference) and may 
transform to any of the three ferroelectric phases. All three phases may have the same free energy, so that there 
is no preferred state. A phase transition from one state to the other with a low activation energy is possible if the 
interface between the two phases are compatible, i.e. if the the deformations and polarization due to the phase 
transitions (i.e. spontaneous strains and polarizations) are compatible across the interface. If the spontaneous 
strains and polarizations are not compatible across the interface, then the compatibility of the kinematic and 
the dielectric displacements have to be satisfied by elastic strains and dielectric displacements, which will give 
rise to internal fields and associated energy states. 

In this work, we will give a brief description of the deformations and polarizations in PZT due to the phase 
transition. We will then introduce the equations of compatibility, which have to be satisfied and will analyze the 
three possible transformations: tetragonal<-*-rhombohedral, tetragonal<->monoclinic, and monoclinic<-»rhombo- 
hedral. 

The micromechanics approach, employed in this work, is that the material may consists of domains with 
different properties (spontaneous deformations and polarizations), which represent the different phases. Within 
the domains, the properties are homogeneous and the material behavior is linear. The interfaces between domains 
are sharp interfaces, where the properties may exhibit stepwise changes. The material may show a nonlinear 
response to applied loads due to moving domain walls or nucleation and growth of new domains. 



2. SPONTANEOUS DEFORMATION AND POLARIZATION 



In the course of a phase transition, the material will be deformed due to changes in the crystal symmety and 
unit cell geometry We will give in brief the deformation gradients due to the phase transitions in question. 
The deformation gradient is the linear transformation which maps reference configuration onto a deformed 
configuration 

x' = Fx 

We distinguish between the deformation due to the phase transition, U, a rigid body rotation, Q, and the general 
deformation gradient, F, (see. Fig. 2). 
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Figure 2. Schematic of a sequence of deformations. U is the deformation gradient due to the phase transition, Q is a 
rigid body rotation. The final deformation gradient, F, which maps a vector x of the reference configuration onto the 
vector x" of deformed configuration is F — QU 

Taking the cubic state as a reference state, the tetragonal, the rhombohedral, and the monoclinic deformation 
gradients due to phase transitions are 

/ a* \ / a r P r (3 r \ I a m S m \ 

U* = a* , U r = \ I3 r a r (3 r , U rn = (3 m 

\ /3* / \ /3 r (3 r a r J \ 7" 1 / 

It should be noted, that these deformation gradients are given in different coordinate systems. While the 
tetragonal and rhombohedral deformation is defined in a coordinate system, which is parallel to principle axis 
of the cubic cell, the monoclinic deformation is given in a coordinate system with principle axis x\ is parallel to 
the pseudo-cubic [110] direction, and x 3 ||[001] pc . To use the deformation gradients on a common basis, either 
the rhombohedral, or the monoclinic deformation gradients have to be rotated by 45° about the [110] pc axis. As 
a result, the monoclinic deformation gradient in pseudo-cubic(ally oriented) coordinates is 

/ \{a m + [5 m ) \{a m -[3 m ) -±5 m \ 
U rn = \{a m -(3 m ) \{a m +[3 m ) ^<5 m , 
V 7 m / 

or the rhombohedral deformation gradient in the rotated coordinate system is 

/ a r + P r V2(3 r 
U r = a r - [3 r 

\ V2(3 r a r 

The polarizations of the tetragonal, rhombohedral and monoclinic state in the global coordinate system parallel 
to the pseudo-cubic axis are 




The parameters of the deformation gradients have been fitted to data according available structural investigations. 
The polarizations can be estimated from deformations and sublattice shifts, but will not be given here. The 
results for the deformations have been plotted in Fig. 3 and details will be published elsewhere. 7 
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Figure 3. Plot of the deformation parameters of the tetragonal, rhombohedral and monoclinic phase of PZT depending 
on the composition across the MPB 

3. COMPATIBILITY OF DOMAINS 

Consider a crystal which volume can be divided into two domains with different deformation gradients (see Fig. 
4). The crystal is considered to be under traction- and field- free boundary condition, so that there are no internal 
fields (stresses and electric fields) in the pre-transformed state, and also no internal fields if the transformation 
is compatible. Given the deformation gradients of the two domains, we want to know, if these two deformations 
can be compatible. If so, what is the normal direction of the interface? In general, there are necessary rigid body 
rotations, Q, to make the two domains compatible (see Fig. 4). 

The deformation gradients of the two domains I and II have to satisfy the following equation of kinematic 
compatibility 8 : 

F 1 - F H = b n 

For a detailed discussion of this equation see for instance Ref. 8. If we allow domain I to be rotated by the 
rotation Q and insert the deformations due to the phase transitions U 1 and U 11 the equation of kinematic 
compatibility becomes 

QU 1 -U 11 = b®n 

b has to be a vector with real components and n is the normal of the interface between the two domains in 
the reference state. If the material is a ferroelectric material, then also the polarizations have to be compatible. 
This means, that the normal component of the polarizations should be continuous across the interface n', i.e. it 
should be compatible after the transformation. Taking into account, that the polarizations are correlated to the 
orientation of the crystal lattice, and due to the deformations the vector n is no longer a normal to the interface, 
the equation of compatibility of the polarization reads 9 : 

(kQP sI -P sII )U u ~ T n = 

where k = ±1, taking into account, that there are usually two opposite polarizations states associated with one 
deformation state. 

Ball and James 10 have found an algorithm, which provides a general solution for the problem of kincmctic 
compatibility, 8 which will be briefly summarized: 

1. Define matrix C = U n ~ V 1 *^ 1 ^ 11 " 1 




Figure 4. Schematic of a crystal divided into two domains with different deformation gradients. To make the two domains 
compatible along an interface with normal n (in reference configuration), or n' (in the deformed configuration) at least 
one of two domains must be rotated by the rotation Q. 

2. If C = I — ► the problem has no solution 

3. If C J, calculate the eigenvalues Ai, A2, A3 of the matrix C. 
Xi > 0, -» Ai < A 2 < A 3 . 

4. A solution exists if Ai < 1, A2 = 1, A3 > 1 

5. Then, there are two solutions: 



1 i ^(l-AQ A!(A 3 -1) 

PV A3 - Ai \ 

K = ±1 

p is choosen so that |n| = 1 and e, are the normalized eigenvectors, corresponding to eigenvalues Ai 

6. The necessary rotation of domain I is obtained by inserting b and n into the compatibility equation: 
Q = (b®n + U 11 ^ U 1 ' 1 

This scheme provides the complete solution for the kinematic compatibility of two domains with arbitrary 
deformations. Therefore, the electromechanical problem which requires both, kinematic and electrical compati- 
bility is over-determined and has in general no solution. It can be shown, that a solution for the coupled problem 
exits only, if the the two domains are symmetry related. 11 I.e., only if the two domains arc different variants 
of the same crystallographic phase, then both the deformations and the polarization can be compatible. This 
is what would be expected and which gives the well known relationships for the non-180 domains in various 
fcrroelectrics. 

On the other hand, this also means that the transition from one ferroelectric phase to another can not be 
fully compatible, no matter which particular transition is considered. For that reason, here we consider only the 
kinematic compatibility. If a compatible interface plane is found, the polarization charges, which appear on an 
electrical incompatible interface, can be calculated (see Sect. 5). 

4. RESULTS 

In order to analyze the possibility of a compatible phase transition, the eigenvalues of the matrices C according to 
the various possible transitions have to be analyzed. A compatible transition is possible, if one of the eigenvalues 
is equal to 1. For the tetragonal<->rhombohedral transition it is found, that this condition is fulfilled, if 



a* = a r - [3 r . 



The tetragonal<->monoclinic transition is compatible, if 

a* = I3 m , 

and the rhombohcdral^->monoclinic transition is compatible if 

/3 m = a r - [3 r . 

The deformation parameters according to these three conditions have been plotted versus the composition 
across the MPB in Figs. 5-7. It is found, that the tetragonal^rhombohedral may be compatible at 
a composition at the very Zr-rich side of the MPB (x ~ 0.54). However, in this compositional region, 
there are no experimental data for the tetragonal phase, so that it may not exist at this Zr-concentrations. 
On the other hand, the tetragonal<-*-monoclinic transition is never compatible according to the accessi- 
ble data and the monoclinic^-s-rhombohedral transition may be compatible at even higher Zr-contents as 
tctragonal<->rhombohedral. 

In general, it can be concluded that one can find an interface between the two domains with comaptible 
deformation, if the out-of-plane deformations in the [110] pc plane are homogeneous. Then, there is a compatible 
interface with a normal vector in the [Il0] pc plane. 

Taking for instance the differences a* = {3 m as a measure of incompatibility for the transition, it can be noted 
that the incompatibility between the tetragonal and the monoclinic phase is almost constant across the MPB, 
but relatively low compared to the other transitions. This may explain, why tetragonal and monoclinic phases 
seem to coexist within the MPB region. 

5. OUTLOOK 

The micromechanical approach to phase transitions in ferroelectrics at the MPB can be extended towards energy 
minimizations for incompatible interfaces. The internal fields arising due to incompatible deformations and 
polarizations can be modeled by distributed interface defects. For polarizations it is obvious, that the charge 
density at the interfaces can be computed from the differences of the polarizations: 

p= (QP sI -P sII )U u ~ T n 
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Figure 5. Plot of the deformation parameters a and a r — f3 r versus the composition of PZT. At the intersection of both 
curves, the tetragonal and the rhombohedral phase could be compatible 
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Figure 6. Plot of the deformation parameters a' and /3 m versus the composition of PZT. As the two curves do not 
intersect, the tetragonal and the monoclinic phase are not compatible 
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Figure 7. Plot of the deformation parameters f3 m and a r — f3 r versus the composition of PZT. At the intersection of 
both curves, the monoclinic and the rhombohedral phase could be compatible 

The analogous model for the kinematic defects are distributed dislocations, expressed by the interface dislocation 
density tensor 12 

In order to calculate the energy contributions of both types of defects and to search a plane with minimal 
energy, it is necessary to know the complete geometry of the domains as well as the the material properties. This 
is far beyond the scope this contribution. 



6. SUMMARY 



The present study presents a micromechanical approach to phase transitions in ferroelectric materials at a 
morphotropic phase boundary. The peculiar property of such kind of materials is that there is the possibility of 
transitions from one ferroelectric phase to another. This transitions have low energy barriers and can then increase 
the degree of freedom to respond on applied loads, if the interfaces between different phases are compatible. The 
compatibility problem has been studied using a well established formalism. It should be noted, that there exist 
in general no solution for the entire electromechanical problem, where both deformations and polarization have 
to be compatible. For that reason, only the kinematic compatibility is considered here. For the particular 
deformation parameters of PZT it is found that the tetragonal^monoclinic transition is very close, but never 
exactly compatible. The tetragonal^rhombohcdral transition has a point of compatibility at x w 0.54 and the 
rhombohedral^monoclinic transition has a point of compatibility at even higher Zr contents. This study helps 
to explain why the tetragonal and the monoclinic phase may coexist in the entire range of the MPB. But all 
phase transitions in question are associated with incompatible phase interfaces which requires to overcome a 
certain critical energy. 

The outlined approach can be used to model the overall behavior of multi- variant and multi-phase crystallites 
with certain, simplified geometrical arrangements of the constituents. 
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